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Centrality dependence of charged antiparticle to particle ratios 
near mid-rapidity in d+Au collisions at y/s NN = 200 GeV 
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The ratios of the yields of charged antiparticles to particles have been obtained for pions, kaons, and 
protons near mid-rapidity for d+Au collisions at ^/s NN = 200 GeV as a function of centrality. The 
reported values represent the ratio of the yields averaged over the rapidity range of 0.1 < y-w < 1.3 
and < Uk,p < 0.8, where positive rapidity is in the deuteron direction, and for transverse momenta 
0.1 < pj' < 1.0 GeV/c and 0.3 < < 1.0 GeV/c. Within the uncertainties, a lack of centrality 



dependence is observed in all three ratios, 
and model calculations. 

PACS numbers: 25.75.-q 



The data are compared to results from other systems 



Experiments at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory aim to 
understand the behavior of strongly interacting matter 
at high temperature and density, testing predictions 
of quantum chromodynamics (QCD). As part of this 
investigation, smaller systems at RHIC energies need 
to be studied in order to aid in the understanding and 
interpretation of results from the more complicated heavy 
ion collisions. In this paper, the ratios of the yields of an- 
tiparticles to particles for primary charged pions, kaons, 
and protons in d+Au collisions at y/s NN = 200 GeV as a 
function of collision centrality were determined using the 
PHOBOS detector during the 2003 run. 

Anti-proton to proton yield ratios near mid-rapidity 
depend largely on the dynamics of baryon-antibaryon 
pair production and baryon number transport in nuclear 
collisions. The rate of pair production can depend on 
the state of the matter created, see references in p). 
Recent comparisons of d+Au and Au+Au data suggest 
that the conditions in Au+Au collisions are very different 
from those observed in d+Au 3, 0, 0, 0- Whether 
these different conditions influence the particle ratios is 
explored by measuring the ratios in d+Au collisions. 

In Au+Au collisions at RHIC energies the ratio of 
(P) I (p) increases from 0.6 to 0.8 as the collision energy 
increases from 130 GeV to 200 GeV, and shows a weak 



dependence on centrality and p T H H H GH El H 
Ibi Il4| . These results imply that baryon-antibaryon pair 
production is larger than baryon number transport and 
yet there is still a finite baryon number transport over 5 
units of rapidity . 

Based on lower energy data, the expectation is that 
the more collisions (v) a participating nucleon suffers, 
the grea ter the baryon number transport to mid-rapidity 
[TH Il6|. This results in the reduction of the (p) / (p) 
ratio. For central Au+Au collisions (12% most central 
events ^3) eacn participating nucleon suffers on average 



5.2 collisions, v = N c °" 2 , where N co u and N part 
are the number of binary collisions and the number 
of participants, respectively. In d+Au collisions, when 
looking in the deuteron hemisphere, v can be defined 
using the number of deuteron participating nuclcons 
(y = ). Over the range of centrality discussed in 

J * part 

this paper, (v) varies from 2 to 8. Hence the range 
of v in this d+Au measurement includes that observed 
in Au+Au collisions. This allows a comparison of the 
relative magnitude of the baryon number transport per 
produced baryon between the two systems. 

Results from Au+Au collisions 0, |n| show that 
the (K~) j (K + ) and (p)/(p) ratios are consistent with 
thermal models. This implies that frequent final state 
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interactions occur. In d+Au collisions little reinteraction 
is expected and therefore the particle ratios reflect the 
initially produced yields. Thus, the influence of baryon 
number transport, baryon production, and final state 
interactions can be investigated by comparing results 
from the two systems. 

The results reported in this paper were obtained using 
the PHOBOS two-arm magnetic spectrometer |ljj. Each 
arm has a total of 16 layers of silicon sensors, providing 
charged-particle tracking both outside and inside the 

2 T field of the PHOBOS magnet. Particles within 
the geometrical acceptance region used in this analysis 
traverse at least 12 of the 16 layers. Three single layer 
silicon pad detectors ( "Ring counters" ) located on either 
side of the interaction point were used to determine 
the multiplicity in the pseudo-rapidity range covering 

3 < M < 5.4. 

Another single layer silicon pad detector ( "OCT" ) sur- 
rounds the interaction region with a cylindrical geometry 
along the longitudinal (z) direction with \z\ < 50 cm, 
corresponding to a pseudo-rapidity coverage of |ry| < 3.2. 
The z position of the vertex is found by maximizing the 
number of OCT hits above a variable threshold. Due 
to the changing angle of incidence, the energy deposited 
per track passing through the OCT detector, and hence 
the low energy cutoff that defines a hit, increases as 
the distance from the vertex increases. The resulting 
vertex position resolution ranges from 0.7 cm to 1.3 cm 
in central and peripheral collisions, respectively. This 
method of vertex reconstruction was found to be the most 
efficient for low-multiplicity events. 

An event-by-event reconstruction of the transverse 
position of the interaction vertex is not possible due to 
the low track multiplicity in d+Au collisions. Instead, 
the average transverse position of the vertex (beam orbit) 
for a given data taking period is used. The beam 
orbit is determined from the intersection points of tracks 
traversing multiple layers of silicon in the spectrometer 
planes that lie outside of the magnetic field and in a two 
layer silicon pad detector covering \r]\ < 1.5 and 25% of 
the azimuthal angle. 

The primary event trigger was provided by two sets of 
10 Cerenkov detectors ("TOs"), which cover the pseudo- 
rapidity range -4.9 < t) < -4.4 (TON) and 3.7 < ?/ < 
4.2 (TOP) for the nominal vertex position (z — 0), where 
positive r\ is defined as the direction of the deuteron. 
These asymmetric positions were chosen to optimize 
the acceptance for primaries, without shadowing the 
Ring counters. A triggered event required a coincidence 
between TON and TOP as well as a time agreement 
that corresponded to an approximate vertex range of 
\z\ < 50 cm. The data set of 30 million triggered events 
requiring only this condition is referred to as "dAVertex" . 
To enhance the sample of peripheral events, a separate 
data set with an additional online trigger condition was 
taken. This trigger required that the occupancy in each 
of the two sets of 16 scintillator paddle counters, which 
cover the pseudo-rapidity range 3 < |?7| < 4.5, be less 



than 50%. This data set of 20 million triggered events is 
referred to as "dAPeriph" . 

Offline event selection cuts were applied to both d+Au 
data sets. To ensure that a non-spurious vertex was re- 
constructed, a more restrictive cut on the time difference 
between the TOs was applied. In addition, a cut was 
applied requiring agreement between the standard and 
TO vertices. To achieve uniformity across all PHOBOS 
trigger configurations, it was required that both sets of 
paddle counters have at least one hit. To reduce vertex- 
position-dependent systematic effects, only events with 
a vertex of \z\ < 8 cm were used. This range was 
chosen to ensure that both particles and antiparticles can 
be tracked and identified in the spectrometer for both 
polarity settings. 

For this analysis, the events were divided into four 
ccntrality classes based on the observed total angle- 
corrected energy deposited (Em ng ) in the Ring counters, 
which is proportional to the number of charged particles 
hitting these detectors. The four centrality classes 
were determined by cuts in Eja ng which correspond to 
a percentage range of the full Em ng distribution not 
biased by trigger or vertex inefficiencies. These cuts 
were determined from a Glauber model calculation using 
HIJING [H and a GEANT 3.21 simulation of the full 
detector. Table[I]shows the percentages that define these 
classes of events. The relative multiplicity of each data 
set to that in the most central 10% bin ((££°™)) is 
also shown to provide a model-independent measure of 
the centrality. Using the HIJING model to relate Em ng 
and N co ii, the average (N co u), impact parameter (b) and 
number of collisions per deuteron participant (is) for each 
trigger condition and centrality bin can be determined. 
In addition, the average trigger and vertex efficiency (e) 
as determined using HIJING is quoted. The additional 
trigger requirements in the dAPeriph data set remove 
more central events, resulting in a reduced efficiency 
and also slightly different average values in the 30-60% 
centrality bin relative to the dAVertex data set. 

The tracking used in this analysis is similar to that 
used in previous Au+ Au analyses 0, El El , and is the 
same as in [2j]. For the low- multiplicity environment 
of d+Au collisions the track seeds can be determined 
without knowledge of the vertex position. This "vertex- 
less" tracking is the major difference between the d+Au 
and Au+Au tracking algorithms. The tracking was 
also extended to include additional spectrometer sensors 
within the field region, as compared to the previous ratio 
analyses 0, El- This results in an increased rapidity 
coverage. A three sigma cut on the track's distance of 
closest approach to the beam orbit (dca< 0.35 cm) was 
applied to reduce the secondary contribution, which is 
largest at low pr. Additional pt cuts of 0.1 GeV/c for 
pions and kaons and 0.3 GeV/c for protons were applied 
to keep the secondary contribution below 1% and 5%, 
respectively. 

Particle identification (PID) was based on the trun- 
cated mean of the specific ionization dE/dx measured 
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FIG. 1: Distribution of average truncated energy loss as a 
function of reconstructed particle momentum. Three clear 
bands can be seen, corresponding to pions, kaons, and 
protons. The solid lines indicate the cuts used for particle 
identification. 
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FIG. 2: Contours of the acceptance of the spectrometer as 
a function of transverse momentum and rapidity for pions, 
kaons, and protons where the raw counts per y and pr bin 
have fallen to 10% of their maximal value. The left plot is for 
particles bending toward the beam pipe and the right plot is 
for particles bending away. The acceptance is averaged over 
the selected vertex range and the accepted azimuthal angle. 

in the silicon spectrometer planes [jj. The PID cuts 
for pions, kaons, and protons are shown in Fig. ^ 
The curved bands are based on the position where the 
(dE/dx) distribution for a given momentum is three RMS 
deviations away from the mean expected value for each 
species. The upper momentum cut for the pions and 
the lower dE/dx cuts for the kaons and protons are 
determined in order to minimize possible contamination 
from other particle species. 

For a given field polarity, particles of both charge 
signs can be reconstructed, but with different kinematic 
acceptances. These different acceptances are outlined in 
Fig. The particles of a given charge bend in the same 
direction and hence have the same acceptance as do the 
particles of the opposite charge in the opposite polarity. 

The acceptance-corrected particle ratios are deter- 
mined for each bending direction simply by the ratios 
of the raw counts per event for antiparticles and par- 
ticles. This procedure assumes that the particle and 
antiparticle acceptance, tracking efficiency and kinematic 



distributions for each bending direction are the same over 
the spectrometer acceptance region for each centrality. 
This was verified by confirming that the opposite po- 
larity field strength and Em ng fractions agree within 
0.2% and 1% respectively. In addition, particle and 
antiparticle average kinematic values {(pt), (pt), (y)) 
agree within 2% for each bending direction. Polarity 
dependent systematic effects that are the same for each 
bending direction are removed by simple averagingof 
the ratios measured for the two bending directions [llj . 
Examples of this include the field strength and centrality 
dependence. Polarity-dependent systematic effects that 
are different for each bending direction, such as vertex 
distributions, must also be taken into account. To correct 
for differences in the beam orbit the data were divided 
into statistically independent subsets. Differences in the 
vertex distribution in the z direction were accounted for 
by applying a z-dependent weight to the raw counts. 

Table [n] gives a summary of the particle and event 
statistics used in this analysis. The raw ratios for 
each centrality bin are determined from the statistically 
weighted average of the ratios over finer subsets of the 
data in order to reduce systematic errors. 

The systematic errors in the raw ratios were deter- 
mined by examining the effects of varying the cuts used 
for event selection, centrality, track selection, and PID 
determination. There were two dominant sources of 
systematic errors. The first originates from the method 
used to determine centrality. Other measures of cen- 
trality lead to a point-to-point systematic error of ±2% 
assigned to the kaon and proton ratios. The second major 
contribution originated from the dependence of the ratios 
on the kinematic acceptance (px and y) over which they 
are measured. An error of ±1% is assigned to the kaon 
and proton ratios. There are additional systematic error 
contributions from dead and hot spectrometer channels, 
spectrometer arm asymmetries, and polarity-dependent 
vertex corrections. For proton ratios they amounted 
to 0.5%, 1%, 1%, respectively. Electron contamination 
was estimated to change the (ir~)/(ir + ) ratio by less 
than 0.1%. All systematic errors on the raw ratios were 
added in quadrature, keeping point-to-point and scale 
systematic errors separate. 

The values of the ratios for detected particles will 
be different from those of particles produced in the 
collision if there is a significant yield of particles from 
secondary interactions and weak decays or a loss of 
particles due to absorption in the detector materials. 
Corrections resulting from additional particle yield were 
described in detail in refs. [3, 1 1 1| . These corrections, 
and their systematic errors, are small because most of 
the unwanted particles can be rejected by tracking in 
the highly segmented silicon detectors which begin only 
10 cm away from the interaction point. For the pion 
and kaon data, the total corrections were estimated to 
be less than 0.5% and 1%, respectively. These values are 
reflected in the final systematic errors of the ratios. 

Corrections to the (p) j (p) ratio are more significant. 
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Trigger 


Centrality 


<e> 


<6> 


(N coll ) 


M 


/ rpnorm. \ 
\^Ring ) 


{*-)/{*+) 


{K~)/{K+) 


(p)/(p) 


Condition 


% 


















dAVertex 


60 - 100 


0.20 


7.4(1.4) 


2.9(1.7) 


2.2(1.3) 


0.14 


0.995 ±0.015 ±0.017 


0.97 ±0.07 ±0.03 


0.84 ±0.04 ±0.04 


dAVertex 


30 — 60 


0.61 


0.6(1.0) 


T A/O A\ 
( .0(3.U J 


^ A/ 1 Q\ 

4.U(l.o J 


A OO 

U.33 


1 AA/1 _l_ A AA*7 _l_ A A1 1 
1.U04 ± 11.00/ ± U.01/ 


A AC _l_ A AO J_ A AO 
O.yo ± 0.03 ± 0.03 


A OA _l_ A AO _J_ A AO 

O.oO ± 0.02 ± 0.03 


dAVertex 


10-30 


0.78 


4.0(1.5) 


12(3.6) 


6.1(1.8) 


0.61 


1.008 ±0.006 ±0.017 


0.97 ±0.02 ±0.03 


0.83 ± 0.02 ± 0.03 


dAVertex 


0-10 


0.84 


3.0(1.4) 


16(4.0) 


8.1(2.0) 


1.00 


1.016 ±0.007 ±0.017 


0.97 ±0.03 ±0.03 


0.86 ± 0.02 ± 0.03 


dAPeriph 


60 - 100 


0.18 


7.4(1.4) 


2.8(1.7) 


2.2(1.3) 


0.14 


0.996 ±0.008 ±0.017 


1.02 ±0.04 ±0.04 


0.86 ± 0.03 ± 0.03 


dAPeriph 


30-60 


0.24 


5.9(1.6) 


6.2(2.7) 


3.7(1.6) 


0.29 


1.014 ± 0.007 ±0.017 


0.97 ±0.03 ±0.04 


0.82 ± 0.02 ± 0.03 



TABLE I: Antiparticle to particle ratios within the acceptance for each centrality bin: (e) is the trigger and vertex efficiency, 
(b) is the average impact parameter, {N co u} is the average number of collisions, {v) is the number of collisions per deuteron 
participant, and (E^™^ is the relative multiplicity as measured by Em ng ; the numbers in parentheses represent the RMS 
of their respective values. The systematic errors in these quantities are 30%, 20%, 15%, 10% in order of increasing centrality. 
The errors on the final ratios represent the statistical and point-to-point systematic errors, respectively. The systematic scale 
errors are not shown. 
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n+ 


K~ 


K+ 


P 
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dAVertex 


60-100 


1004052 


2787 


13879 


163 


537 


310 


556 


1120480 


15320 


3087 


524 


157 


543 


421 


dAVertex 


30-60 


2863440 


14963 


74941 


824 


2732 


1604 


3494 


3078602 


81368 


16164 


2930 


1011 


2998 


2282 


dAVertex 


10-30 


2444774 


19663 


95755 


1113 


3660 


2198 


4558 


2612126 


103027 


20641 


3780 


1237 


3824 


2903 


dAVertex 


0-10 


1291300 


13628 


64357 


803 


2497 


1549 


3077 


1378206 


69664 


14303 


2675 


864 


2649 


1958 


dAPeriph 


60-100 


3021106 


8255 


40530 


455 


1416 


866 


1841 


3222936 


42931 


8797 


1509 


461 


1564 


1085 


dAPeriph 


30-60 


3053872 


13666 


65799 


737 


2353 


1461 


3182 


3311696 


71542 


14435 


2449 


821 


2707 


2009 



TABLE II: Summary of particle and event counts. 
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FIG. 3: Particle ratios as a function of centrality for each 
species. The open symbols are from d±Au collisions, the 
filled symbols are from central Au±Au collisions |ll| . both at 
y/s NN = 200 GeV. The brackets represent the point-to-point 
systematic error. The systematic scale errors are not shown. 



These were determined as in |ll| with several important 
differences. The dominant change to this ratio is due to 
the absorption correction, which arises from an asymme- 
try in the loss of anti-protons versus protons interacting 
in the beam-pipe and planes of the spectrometer. The 
correction value of 3.5% ± 1.4% (syst.) was determined 
using GEANT with two hadronic interaction packages, 
Gheisha and Fluka. The systematic error represents half 



of the difference between the two interaction packages 
and is a scale error for all of the proton ratios. Using 
HIJING events, the secondary correction for the produc- 
tion of protons in the beampipe and detector materials 
was found to be 1.6±0.3% (stat.). The correction is larger 
than in the Au±Au analysis [llj due to the different 
acceptance and the reduced ability to define the event 
vertex. This correction is dependent on the HIJING p/it 
ratio matching the true p/it ratio. It was found that this 
assumption is correct at the 10% level, which resulted in a 
systematic scale error of ±0.2% in the correction. A feed- 
down correction of —0.5% accounts for the difference in 
the number of hyperons (primarily A and A) decaying 
to protons versus antiprotons. The reduced ability 
to determine an event-by-event vertex requires that a 
momentum dependent correction be employed. The anti- 
hyperon to anti-proton ratio (A/p) was estimated to be 
0.6 ± 0.3 (syst.), derived from the ave rage of two 
extreme limits, the published Au±Au ratio [20l l2lT | and 
the HIJING d±Au ratio, where the model is known to 
underestimate strangeness production. The systematic 
error reflects the difference between the two values. The 
resulting point-to-point systematic error in the correction 
for the two most central bins is ±1% and for the two 
peripheral bins is ±0.6%. 

Table [I] shows the final antiparticle to particle ra- 
tios within our acceptance for the different centrality 
bins. The systematic scale error on the ratios has 
been separated from the point-to-point errors, shown in 
Table [J allowing for a more precise determination of the 
centrality dependence of the ratios. The systematic scale 
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FIG. 4: Antiproton to proton ratio comparison with models. 
The squares are the d+Au data at ^/s NN — 200 GeV, with 
brackets representing the point-to-point systematic error. 
The lines are fits to the model predictions of (p)/{p) within 
the acceptance. The statistical error in the models is less than 
2%. 

errors for pions, kaons and protons are ±0.008, ±0.02 
and ±0.02, respectively. For the proton ratio a large 
contribution to this error comes from the uncertainty in 
determining the absorption correction. 

Fig. |31 shows the final ratios for each species as a 
function of centrality. For clarity of presentation, the 
two measurements where the two data sets overlap are 
weighted together statistically. Within the statistical 
and systematic errors all the particle ratios appear to 
be independent of centrality in d±Au collisions. The 
12% central Au±Au ratios 0] are also shown in Fig. 
|3J The value of (v) « 5.2 for central Au±Au collisions 
is determined from a Glauber model calculation using 
HIJING. The pion and kaon ratios agree between the two 
systems. This suggests that any final state interactions in 
Au±Au do not modify the ratio of the initially produced 
meson yields. In contrast, the Au±Au proton ratio is 
significantly lower than the ratios in all of the centrality 
bins in the d±Au collisions. Additional data from 



Au+Au collisions BEHU 

suggest that for less central 
Au+Au collisions (lower values of v) the (p)/(p) ratio in 
our acceptance is similar to that found in the d±Au data. 

Assuming that the total proton yields are the sum 
of a transported component and a produced compo- 
nent which is equal to the antiproton yield, the value 
(p)/(P) ~ 1 is a measure of the relative fraction of 
transported protons to produced protons. A comparison 
of the central d±Au and Au±Au results .11] show that 
the relative fraction of transported protons in a central 
d±Au collision is half that observed in a central Au±Au 
collision, despite the larger value of v in the central 
d±Au collisions. This may be evidence of collective 
behavior that affects baryons in Au±Au collisions and 
is not present in d±Au collisions. 

Fig. 0] compares the (p)/(p) ratio as a function of 
centrality with the same ratio from HIJING [18J , RQMD 
and AMPT jHHI. The model outputs were passed 
through a simulation of the detector and the same trigger, 
event, and particle selection biases as used in the data 
analysis were applied. These models agree with the 
expectation of increased baryon transport with increasing 
v, which results in a decreasing ratio. The ratios averaged 
over different centralities from the models and the data 
are roughly consistent. However, suppression of the 
ratio with increasing centrality seen in the models is not 
observed in the data. 

In conclusion, the data shown in this paper provide 
the first information on the baryon transport in d±Au 
collisions at the full RHIC energy. These ratios provide 
constraints on current and future models dealing with 
baryon production and transport and thus set a baseline 
from which to further investigate Au±Au and other 
systems. 
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